Screening and finite size corrections to the octupole and Schiff moments 
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Parity (P) and time reversal (T) violating nuclear forces create P, T-odd moments in expansion 
of the nuclear electrostatic potential. We derive expression for the nuclear electric octupole field 
which includes the electron screening correction (similar to the screening term in the Schiff moment). 
Then we calculate the Z 2 a 2 corrections to the Schiff moment which appear due to the finite nuclear 
size. Such corrections are important in heavy atoms with nuclear charge Z > 50. The Schiff 
and octupole moments induce atomic electric dipole moments (EDM) and P, T-odd interactions in 
molecules which are measured in numerous experiments to test CP-violation theories. 
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The measurements of EDM induced by parity- and 
time- violating forces provide crucial tests of modern uni- 
fication theories. The best limits on P, T-odd nuclear 
forces are obtained from the measurements of 199 Hg 
atomic EDM Nuclear EDM can not induce atomic 
EDM due to complete electron screening of nuclear EDM 
(Schiff theorem 03 )• The first P, T-odd terms that 
survive the screening and produce atomic EDM are the 
Schiff and octupole moments [5]. In this paper we con- 
sider the partial screening of the nuclear electric octupole 
moment and the finite nuclear size corrections to the 
Schiff moment. 

The nuclear electrostatic potential with electron 
screening taken into account can be written in the fol- 
lowing form (see e.g. @ for the derivation): 
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where J p(r)d 3 r — 1, and d = Ze(r) = Ze J p(r)rd 3 r is 
the nuclear EDM. The second term cancels the dipole 
long-range electric field in the multiple expansion of 
<^(R). We expand the Coulomb potential in the terms 
of Legendre polynomials 
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where r< and r> are min[r, P] and max[r, R] respectively, 
6 is the angle between vectors r and R. The Legen- 
dre polynomials can be written in the following form: 



Pi (cos 0) = nRi/(rR), P 2 (cos#) = P ( P, -q tJ /(2r 2 R 2 ) 
where = ir^j — r 2 6ij is the quadrupole moment ten- 
sor (summation over repeating indexes is assumed) . The 
P, T-odd part of the potential ([1]) originates from the odd 
harmonics I of the first term and even harmonics of the 
second term. We start our consideration of P, T-odd part 
ip {1) {R) of ([I]) with I = 1 in the first term and I = 0,2 



in the second term. The dipole part of (p^-'(R) corre- 
sponds to the Schiff moment field. The third harmonic 
/ = 3 in the first term of |T]) gives the octupole field that 
has been considered in Q- We will add I = 3 later. As 
we will show below, account of I = 2 in the second term 
of |T]) gives the screening of the octupole field. We can 
present the P, T-odd part of the potential ((T|) in the form 
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Note that for R —> oo the first and third terms of Eq. §3§ 
cancel each other. Therefore, for their sum we can use 

Jcf = Jo°° - Ir = - Ir and Present </>W as 
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Last term in the above equation can be presented as 
{rj)Rj 5 (r k ) RjRjRk _ 5 
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In Eq. ([5]) the tensor in the brackets {...} vanishes after 
convolution with the symmetric tensor qij in Eq. 
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Introducing Oij k into Eq. (j4]) we obtain the equation for 
the P, T-odd part of the electrostatic nuclear potential 



5 Ze 



p(r)d 3 r 



The last term in the above expression originates from the 
second (screening) term in Eq. ([TJ . It gives the screening 
for the octupole field. The octupole appears due to the 
third harmonic I = 3 in the Coulomb potential expansion 
and for R > Rn is given by the equation: 
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In the above equation the nuclear octupole moment ten- 
sor is 
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Taking the screening term into account we can present 
the screened nuclear octupole moment tensor in the fol- 
lowing form: 
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{(n)q 3 k + (rk)qi 3 + (rj)qu} ■ (10) 



We see that the octupole screening is expressed in terms 
of the nuclear electric dipole moment (d = Ze(r)) and 
the nuclear quadrupole moment operator The par- 
tially screened octupole potential is given by 
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The inner part of the nucleus does not give significant 
contribution to the electron matrix elements of the nu- 
clear octupole field. The situatution is different for the 
Schiff moment field which was considered in Refs. 
To make the picture for the electrostatic T, P-odd nuclear 
potential complete, we will present a brief derivation for 
tthe Schiff moment field including the finite nuclear size 
corrections. The Schiff moment field is given by the first 
term in Eq. 0), 



(ls \K) = ZeR, 
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p(r)d 3 r. 
(12) 

We see that v? (15 '(R) = if R > Rn (nuclear radius) 
since p(R) = in that region, i.e. this potential is local- 
ized inside the nucleus. 

The electric field of the nuclear Schiff moment polar- 
izes the atom and produces atomic EDM. All electron 
orbitals for I > 1 are extremely small inside the nucleus. 
Therefore, we should consider only the matrix elements 



between s and p Dirac orbitals. We will use the following 
notations for the electron wavcf unctions: 
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where Clji m is a spherical spinor, n = R/i?, f(R) and 
g{R) are the radial functions. Using (<x ■ n) 2 = 1 we can 
write the electron transition density as 



p sp (R) = = ^pU sp (R) 



(14) 



U sp (R) =f s (R)f p (R) + g s (R)g P {R) = V b k R k (15) 
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The expansion coefficients b k can be calculated analyti- 
cally [8|; the summation is carried over odd powers of k. 
Using Eqs. (|12I14I) we can find the matrix elements of 
the electron-nucleus interaction, 
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where (s\n\p) = J Qt n n p # sin Qd6, (r n ) = J p(r)r n d 3 r. 
Note, that all vector values (rr n ) are due to P, T-odd 
correction bp to the nuclear charge density p a , while (r n ) 
are the usual P, T-even moments of the charge density 
starting from the mean-square radius (r 2 ) = r 2 for k = 1. 

In the limit of the point-like nucleus the Schiff moment 
potential and its matrix element are given by [5| : 

V?s(R) =4vrS- V5(R) (17) 

(a|-e Vs |p)=47reS-(V^ P )fl=o. (18) 

For the solutions of the Dirac equation the product 
(V-0jV'p)fl-yo is infinite for a point-like nucleus. There- 
fore, we need a finite-size Schiff moment potential. We 
have shown Q that this potential increases linearly in- 
side the nucleus and vanishes at the nuclear surface. We 
suggested an approximate expression for such potential 
which is convenient for the calculations of atomic EDM: 



, 3S • R 

<^ S (R) = —n(R), 



(19) 



where B = J n(R)R A dR « R%/5, Rn is the nuclear 
radius and n(R) is a smooth function equal to 1 for 
R < Rn ~ S and for R > Rn + 5; n(R) can be taken as 
proportional to the nuclear density po (note that we can 
choose any normalization of n(r) since the normalization 
constant cancels out in the ratio n/B, see Eq. (fT9"]0 . 
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TABLE I: Schiff moment S and the ratio of the Za correction AS = 5" - S to S for 205 T1 and 199 Hg nuclei. So and A So are 
the bare values, without the core polarization correction produced by the strong residual nuclear forces. Stot and AStot are the 
results including the core polarization effects. In brackets there are values of L'/S ( where L' — L — S) which are the results 
of the calculations from To obtain the final values of S and S' — S + AS one should sum the contributions of the three 
interaction constants go, g\ and The units of S are e • /m 3 . 







So 


AS/So (L'/So) 


Stot 


AStot/ Stot (L'tot/ Stot) 




9a 


-0.085 


0.05 (-0.1) 


-0.006 


0.09 (-0.05) 


199 Hg 


9i 


-0.085 


0.05 (-0.1) 


-0.036 


(-0.15) 




92 


0.17 


0.05 (-0.1) 


0.019 


0.06 (-0.08) 




90 


-0.075 


0.05 (-0.09) 


-0.014 


-0.03 (-0.18) 




9i 


-0.028 


0.23 (-0.39) 


-0.082 


-0.03 (-0.18) 




U-2 


0.237 


0.06 (-0.08) 


-0.007 


-0.34 (-0.51) 



We now set the matrix elements (fTpjl of the true nu- 
clear T, P-odd potential (|12p to be equal to the matrix 
elements of the effective potential (|T9| which are given 
by 

(s| - ecf>(R)\p) = 15e(«|n|p> ■ -g- / U sp R 3 n(R)dR 

N JO 



6 3 /6i « -0.5Z 2 a 2 /R 2 N . This gives 
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where we have made approximation 
/ n(R)R k dR w i?^ +1 /(fc + 4). Equating ^T6j) and 
(l20l we obtain 



S' = 
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Note that S' in Eq. (|21l) differs from the local dipole mo- 
ment L defined in Ref . [8j and calculated in [9j . The local 
dipole moment L does not contain the sum in the de- 
nominator in Eq. (|21j) and corresponds to the <5-function 
form (similar to Eq. (fTT)) ) of the effective Schiff moment 
potential. 

In light atoms (Za <C 1) it is sufficient to keep 61 only. 
This gives us the well-known expression for the Schiff 
moment 



This equation allows one to calculate the Za corrections 

(20) to the Schiff moment. In our work Q such calculations 
were performed for the local dipole moment L. Corre- 
sponding expression for the local dipole moment L d, Q 
does not contain the factor (1 — 5Z 2 ev 2 /14) in the de- 
nominator (see Ecl (|!Z3"]) for 5'). All recent atomic EDM 
calculations |l]] - tl3| have been performed using the Schiff 
moment potential in the form (|19[) . Therefore, we calcu- 
late the values of S' for two nuclei of experimental inter- 

(21) est, 199 Hg (with valence neutron) and 205 T1 (with valence 
(r.-)(q ■ r k ~ 1 ) proton), using calculations of L in Ref. [t| . 

In [9( we used a finite-range P, T- violating nucleon- 
nucleon interaction of the form 
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{[So^a • T b + 52 (r a • r b - 3t^t 6 2 
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The first correction is given by the ratio 63/61. 
This ratio is different for matrix elements s - 
P1/2 (63/61 = -(Zlh)Z 2 a 2 jR 2 N ) and s - p 3/2 
(63/61 = -(9/20)Z 2 a 2 /i? 2 r ). However, with the 10% 
accuracy we can use the average of these two values 



where m p is the proton mass and r au — |r„— r&|. The core 
polarization corrections produced by the strong residual 
nuclear forces have been calculated using the RPA tech- 
nique. 

Results of the calculations of the Schiff moment and 
Za corrections are presented in Table HI As one can see, 
in most cases use of S' (instead of L) leads to smaller 
values of the Za corrections. For S' typical values of the 
Za corrections are about 5-10%. Larger Za corrections 
(up to 34% for S' and 51% for L) appear in the cases 
where the main contributions are suppressed. 
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